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Abstract 

We extend the Rome-Southampton regularization independent momentum- 
subtraction renormalization scheme(RI/MOM) for bilinear operators to one 
with a nonexceptional, symmetric subtraction point. Two-point Green's func- 
tions with the insertion of quark bilinear operators are computed with scalar, 
pseudoscalar, vector, axial-vector and tensor operators at one-loop order in 
perturbative QCD. We call this new scheme RI/SMOM, where the S stands for 
"symmetric". Conversion factors are derived, which connect the RI/SMOM 
scheme and the MS scheme and can be used to convert results obtained in 
lattice calculations into the MS scheme. Such a symmetric subtraction point 
involves nonexceptional momenta implying a lattice calculation with substan- 
tially suppressed contamination from infrared effects. Further, we find that 
the size of the one-loop corrections for these infrared improved kinematics is 
substantially decreased in the case of the pseudoscalar and scalar operator, 
suggesting a much better behaved perturbative series. Therefore it should al- 
low us to reduce the error in the determination of the quark mass appreciably. 



1 Introduction 



Lattice simulations in quantum chromodynamics (QCD) allow for ab initio nonpertur- 
bative determinations of operator matrix elements and physical quantities such as quark 
masses and the strong coupling constant. One starts with a direct computation of the 
bare quantities with the lattice spacing acting as the ultraviolet cutoff in some particular 
discretization of QCD. Providing that the lattice spacing is sufficiently small, it is in 
principle possible to obtain the corresponding renormalized quantities using perturbation 
theory. However, the coefficients in lattice perturbation theory frequently prove to be 
large and for this reason techniques using nonperturbative renormalization (NPR) have 
been developed and are being successfully implemented. With these techniques lattice 
perturbation theory is avoided entirely, and one obtains renormalized quantities in some 
appropriate renormalization scheme such as the regularization independent momentum- 
subtraction (RI/MOM) scheme pp. 

On the other hand perturbative calculations in continuum QCD are conventionally and 
conveniently performed using dimensional regularization [2] and the MS renormalization 
scheme [31 H] which is not directly amenable to the NPR procedure. The continuum 
perturbation theory is therefore used to match the quantities computed in the RI/MOM 
and MS schemes. For example the computation of the mass conversion factor Cffl MOM , 
which converts a quark mass renormalized in the RI/MOM scheme into the MS scheme or 
the conversion factor (7^/ MOM , which performs the corresponding conversion of the quark 
fields, are both known up to three-loop order in perturbative QCD [HEHH]. Another 
scheme, which is useful in lattice simulations is the RI'/MOM scheme in which these 
conversion factors are also known up to three- loop order [6j[7] . A more detailed definition 
of these schemes will be discussed in Section EJ The conversion factors C m and C q in both 
schemes can be obtained through the evaluation of self-energy diagrams. Not only quark 
masses, but also the strong coupling constant a s has been studied in MOM schemes [8HT3]. 

With regard to the vertex diagrams one has many choices of defining the subtrac- 
tion point at which the renormalization constants are fixed through different momentum 
configurations. In this paper we determine the one-loop matching coefficients for a gener- 
alization of the RI/MOM scheme in which there are no channels with exceptional momenta 
and which we propose to use in our numerical simulations. Because the kinematic con- 
figuration in this scheme is symmetrical in the three channels, we call it the RI/SMOM 
scheme. In the following we define the symmetric and asymmetric Minkowski momentum 
configurations by 

• symmetric or nonexceptional momentum configuration: 

Pi = p\ = Q 2 = ~/A /i 2 > 0, q = Pi- p 2 , 

• asymmetric or exceptional momentum configuration: 

P 2 i=pI = -V 2 , V 2 > 0, Pi=P2, q = Q, 
where the momentum flow is shown diagrammatically in Fig. [TJ 

In Ref. [H] quark masses were determined through lattice simulations using nonper- 
turbative renormalization [1] in the RI/MOM scheme and subsequently converted to the 
MS scheme. In order to renormalize the bare quark masses in the lattice simulation, the 
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Figure 1: Momentum flow of a generic diagram required for the renormalization procedure 
with nonexceptional momenta. The gray bubble stands for an operator insertion and 
higher order corrections. 

renormalization constants need to be computed on the lattice. In regularization and renor- 
malization schemes which preserve flavor and chiral symmetries in the limit of vanishing 
quark masses, the perturbative renormalization constants of the axial-vector and vector 
operators as well as the ones for the pseudoscalar and scalar operators need to be equal. 
In the standard RI/MOM and RI'/MOM schemes the normalization conditions for quark 
bilinear operators are imposed on Green's functions with the operator inserted between 
equal incoming and outgoing momenta say p, and — p 2 = fi 2 is the renormalization scale. 
The momentum q inserted at the operator is therefore so that there is an exceptional 
channel, i.e. one in which the square of the momentum is much smaller than the typical 
large scale (/i 2 ). For the asymmetric subtraction point effects of chiral symmetry breaking 
vanish only slowly like 1/p 2 for large external momenta p 2 . In Ref. [15J it was proposed 
instead to use a similar renormalization procedure but with the incoming and outgoing 
quarks having different momenta, p\ and P2 respectively, with p\ = p\ = (pi — p^) 2 = p 2 . 
There are now no exceptional channels and we explain below that this decreases chiral 
symmetry breaking and other unwanted infrared effects. The choice of such a symmetric 
subtraction point is very convenient, the renormalized quantities depend also only on a 
single scale p 2 . When the renormalization constants of quark bilinear operators are fixed 
at a symmetric subtraction point (chosen to have nonexceptional kinematics) chiral sym- 
metry breaking and other unwanted infrared effects are better behaved and vanish with 
larger asymptotic powers of the order 1/p 6 . This behavior has been derived in Ref. [15] 
as a consequence of Weinberg's theorem [16] and demonstrated by explicitly computing 
the renormalization constants on the lattice. Hence these RI/SMOM kinematics suppress 
infrared effects much more strongly than the usual exceptional configuration for large ex- 
ternal momenta. The symmetric momentum configuration is thus much more favorable. 
However, in order to be able to use it to evaluate the matrix elements of quark bilinear 
operators and the quark mass, the matching factors need to be determined perturbatively 
for this new, symmetric choice of momenta. A nonperturbative test of the RI/SMOM 
scheme for the quark mass renormalization can be found in Ref. [T7] . 

Another drawback in the case of the exceptional momenta is that the perturbative 
expansion of the usual conversion factor C* I/MOM shows poor convergence and makes a 
significant contribution to the systematic uncertainty in the quark masses obtained from 
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the lattice studies. In fact, in Ref. [H] the error (« 11%) in the quark masses arising from 
the truncation of the perturbative series in the matching factor amounts to around 60% of 
the total error. Therefore determining the conversion factor for a symmetric momentum 
configuration will also allow us to see if the convergence will be better behaved. If it is 
better behaved, then the symmetric configuration would be preferred for both of these 
reasons. Motivated by these considerations we study in this work the renormalization 
of quark bilinear nonsinglet operators of the form O = uTd for a symmetric subtraction 
point, where V represents a Dirac matrix and u and d are fermion quark fields. 

Even with the use of the symmetric, nonexceptional kinematics, the renormalization 
prescription is not unique and the chiral Ward-Takahashi identities can be satisfied us- 
ing a variety of procedures. In the following sections we study a specific scheme which 
we consider to be convenient and practicable for the nonperturbative renormalization of 
lattice quark bilinear operators. In order to preserve the Ward-Takahashi identity, the 
definitions of the vertex and wave function renormalizations are related as we explain in 
the following Section. 

The outline of this paper is as follows: In Section [2] we define our notation and con- 
ventions and introduce the framework required for performing renormalization of the 
quark bilinear operators with a symmetric subtraction point. Subsequently we present in 
Section [3] two methods for the extraction of the conversion factor C m in the RI/SMOM 
scheme, apply the concepts of Section [2] to calculate the vector, axial- vector, pseudoscalar, 
scalar and tensor operators between two off-shell quark states at one-loop order in pertur- 
bative QCD for the nonexceptional momentum configuration and determine the matching 
factors. Finally we close with a brief summary and our conclusions in Section HI Even 
with the symmetric nonexceptional kinematics the choice of renormalization conditions 
is not unique. In Appendix |A] we therefore present the one-loop perturbative results in a 
form which can be used to calculate the conversion factors from a general scheme with a 
symmetric subtraction point to the MS scheme. For illustration we study one alternative 
scheme called the RI/SMOM^ scheme, in which the vertex renormalization condition 
is the same as in the RI/MOM scheme, but with nonexceptional kinematics and with a 
different wave function renormalization. We also provide the results for the conversion 
factors and, in Appendix (Bj the corresponding two-loop anomalous dimensions. 

2 Concepts and framework of the RI/SMOM scheme 

We will begin with a bare, continuum theory of QCD which has been regulated using a 
scheme which guarantees that Green's functions involving the quark field and quark field 
bilinears obey the usual chiral and flavor symmetries of QCD. Dimensional regularization 
is an example of such a scheme. 

Let us consider the nonamputated Green's function Gq of an operator O computed 
between two external off-shell quark lines in a fixed gauge. The corresponding diagrams 
up to one-loop order in perturbative QCD are shown in Fig. [2J 
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(a) (b) (c) (d) 

Figure 2: Diagrams contributing to the nonamputated Green's function up to one-loop 
order in perturbative QCD. The black box indicates the inserted operator. Spiral lines 
denote gluons and solid lines fermions. 

The amputated Green's function is defined by 

A 6 = S- 1 ( P2 )G 6 S' 1 ( Pl ) (1) 
where S(p) is given by the quark propagator 

" = /^<T[* W *(0)]> = ,_„, + ' fe _ E(p) , (2) 

where contains the higher order corrections and can, in perturbation theory, be 

decomposed into its Lorentz structure: E(p) =p , T l y(p 2 ) + rriEs(p 2 ). The lowest order and 
one-loop diagrams contributing to E(p) are shown in Fig. [3] 



p 




(a) (b) 

Figure 3: Propagator- type diagrams up to one-loop order in QCD 

In the following we will consider quark bilinear operators O = uTd with scalar (r = 1), 
pseudoscalar (r = 275), vector (r = 7^), axial- vector (r = 7^75) and tensor (r = 
a fiu _ t^yH^v-^ kernels. We will distinguish between bare and renormalized quantities 
by assigning the index b to a bare quantity and the index r to a renormalized one. In the 
case of renormalized quantities an additional quantifier specifying the scheme is attached. 
Renormalized and bare quantities are related through the renormalization constants Z: 

* R = ZV 2 * B , m R = Z m m B , 6 R = Z 6 6 B . (3) 

The renormalization constants of the scalar (O = S), pseudoscalar (O = P), vector 
(O — V), axial- vector (O = A) and tensor (O = T) operator will be denoted as Z s , Z P , 
Zy, Za and Zt, respectively. In the RI/MOM scheme the renormalization conditions 
which fix the renormalization constants Z m and Z q are given by 

'^ aS *W ] = _i ( 4 ) 



lim 



1 



m R ^o 12m R 



T*[Sr l (p)) 



and lim — Tr 

m R ^0 48 



4 



where the symbol "Tr" denotes the trace over color and spins. The second equation 
determines ^/ MOM and subsequently the first one can be used to extract Z^ MOM . Now 
in the RI'/MOM scheme the second condition of Eqs.fllJ) is replaced by 
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1. 



(5) 



The quark propagator in the RI'/MOM scheme is fixed to its lowest order value at the 



point p 



-fj, , where p is the squared, external, Minkowski momentum and /x is the 



renormalization scale. 

The propagator and vertex diagrams (Fig. [2} for the vector and axial-vector operators 
are related through the vector Ward-Takahashi identity for degenerate masses m u = m d = 
m 

Q»K,b(Pi,P2) = S B 1 {p 2 ) - S^&i) (6) 
and the axial-vector Ward-Takahashi identity 

ilhS B 1 {pi) - S B 1 (p 2 )i^ (7) 
P2- The renormalized and bare amputated Green's 



- «9/iA£ )B (pi,p 2 ) = 2m B Ap !B (p 1 ,p 2 ) 



with the momentum transfer q 
functions are connected by 



Pi 



S R {p) = Z q S B (p), A dR {p 1 ,p 2 ) 



In the following we want to renormalize the quark bilinear operators using a symmetric 
subtraction point. For functions /, which are restricted to the symmetric momentum 



configuration we use the shorthand f(j?i,p%, q 



= f(pi,P2,q 



| sym 



and for the 



asymmetric subtraction point we introduce the abbreviation f(pi,p%, q )| g=0)P 2 = _ A( 2 = p2 = 

asym- 

We perform the quark mass and wave function renormalization by imposing on the 
two-point function S(p) the condition of Eq.© and 



lim 

m R -¥0 I2mn 



Tr 



Sn'ip) 



-ft- 



2 2 



Q^A,R(,Pl,P2)l5 



(9) 



sym 



The second term in the curly brackets on the left-hand side of Eq. (JHJ) starts at 0(a s ) 
and is absent in the RI/MOM and RI'/MOM schemes. This term is needed to maintain 
the Ward-Takahashi identities for renormalized quantities, as we will see below. For the 
vector and axial-vector quark bilinear operators we impose the conditions 



1 



lim 

m R ->o \2q 2 



Tr 



%K,r{p^p*) i 



i. 



sym 



lim 

7iH-K) I2q 2 



Tr 



^ A a,r(Pi,P2)75 



1. (10) 



sym 



The projectors for the amputated Green's functions in Eqs. fllOl) are different from those 
used in the RI/MOM scheme(see Tabled]). Using instead these original RI/MOM pro- 
jectors leads to a different wave function renormalization and will be discussed in Ap- 
pendix |Al For the pseudoscalar and scalar amputated Green's functions we use the 
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renormalization conditions 



lim — Tr [A PR (p 1 ,p 2 )j 5 ] 



1. 



sym 



lim n 7^ Tr 1 A s,r(p^P2) 1 } 



sym 



and for the tensor operator the condition 

1 



lim 

m R ^0 144 



-Tr 



1. 



:n) 



(12) 



sym 



Note that all of the renormalization schemes being considered in this paper are mass- 
independent. Thus, each condition is imposed at fixed external momentum and vanishing 
quark mass. The renormalization conditions of the RI/MOM and RI/SMOM schemes are 
summarized in Table [TJ 



RI/MOM 



lim i Tr 

m R ^0 48 

lim ^Tr 

m R -^0 4S 



A y,R(Pl^2)7M 



= "1, 



m R ^0 4 « L " " 

m hm o ±Tr{AsAPuP2)l] 

v 1 rrLrrr-1 /- \ A 



lim ^— Tr^Vp)] 



= 1. 

asym 



lim ^Tr 

m R ^0 4S 



A' 



hR (Pl,P2)75lti 



= 1, 

asym 



lim ^TrLSo 1 ^)^! , " = -1, 



lim n lk Tr [Vfi(Pl^2)75 



= 1, 

asym 



= 1. 

asym 



RI/SMOM 



lim 1 



- | Tr 

p2 — _^2 I 



%A-A, R {PUP2)1Z 



= 1, 



n 12V Tr k A v,ii(Pi)P2) l/l = 1, hm ^Tr [g^A^ ifl (pi,p 2 )75 4 

m R ^0 lz( < L ^ ' J sym m R ^0 ±zy L ^ ' 

lim n T2 Tt l A sAPi jPs) 1 ] = !> lim n lk Tr [ A P,fi(Pi'P2)7; 



= 1, 

sym 



= 1. 

sym 



Table 1: The renormalization conditions for the RI/MOM and RI/SMOM schemes. 

In the remainder of this Section we will show that if the normalization conditions in 
Eqs.Q and (}9l)- ({TTI) of this RI/SMOM scheme are imposed on the quark bilinear oper- 
ators, the Ward-Takahashi identities of Eqs. flS} and (J7J) are also obeyed for the resulting 
renormalized quantities and the properties Zy = 1 = Zp = 1/Z m and Z$ = Zp are 
preserved, as they are in the MS, RI'/MOM and RI/MOM schemes (see e.g. Ref. [151115] ). 
Some of these properties hold nonperturbatively while the others are proven only in the 
perturbation theory as we will see below. 



Let us start by considering the object y^p-Tr^A^ d\\ S ym and insert the vector Ward- 
Takahashi identity of Eq. (jH]) : 



12q 



sym 



^) rf -Ti[^) 4\ sy z-w ns * 1{q) *C ■ (13) 



12q 



Expressing bare quantities in terms of renormalized ones using Eq. (j8j and imposing the 



condition in Eq.(|5]) and the one on the left in Eq.l flOl) leads to Z 



larly one obtains Z v = Z 



RI/SMOM 



RI/SMOM 
V 



1. Simi- 



by inserting Eq.© into I ^Tr[g At A^ B 7 5 



| symi 
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combining it with Eqs. (fT3|) and imposing the conditions of Eq. (flO|) for the renormalized 
quantities in the massless limit. Note that the above derivation of ^ I / SMOM — ^m/smom 
is independent of the choice of the renormalization point fi. This is in contrast to the 
RI/MOM scheme for which the Ward-Takahashi identity for the axial current only holds 
at large /x 2 . The renormalized vector current satisfies the Ward-Takahashi identity in 
both the RI/MOM and the RI/SMOM schemes even in the low energy region. However, 
the relation Z& = Zy = 1 implies that the axial vertex function given in Eq. ( flUj) remains 
exactly equal to one in the limit of vanishing quark mass even when evaluated in the 
infrared region of QCD where large vacuum chiral symmetry breaking might have been 
expected to introduce large asymmetries between such vector and axial-vector correlation 
functions. 

From ^Vsmom _ ^ _ ^ri/smom ^ f }} ows that the renormalization constant ^/smom 
can be extracted from Eqs. ffTOl) . However, since Eq.([5]) determines Z q in both the RI/SMOM 
and RI'/MOM schemes, zfV SMOM = zf v / MOM , whose value is known up to order a 3 s in 
Ref. [SIC]- Nevertheless, in Section \'d.l\ we will renormalize the vector and axial- vector 
operators for the symmetric momentum configuration in the RI/SMOM scheme using 
the conditions in Eqs. ffTtJj) in order to demonstrate that the value of ^m/smom obtained 
from Eq. ffTO]) is in fact equal to the value for ^m'/mom obtained from Eq.([5]) by explicit 
calculation up to one-loop order. 

From the axial Ward-Takahashi identity it follows that the renormalization constant 
for the pseudoscalar operator Zp l ^ SMOM ^ e mass renormalization constant ^m/smom 
are related. If one multiplies Eq.dTJ) by (—275), takes the trace of both sides over spin and 
color and restricts it to the symmetric momentum configuration, one obtains 



— T:[g M A^ B 7 5 ] 



= 2m B — Tr[A P;B 7 5 ] 

sym lZ% 



sym 6 



1 Tr^)- 1 ! 2 . (14) 



Taking the zero-mass limit, expressing again the bare equation with the help of Eqs.(j3]) 
and OH]) in terms of renormalized quantities and imposing the conditions in Eqs.flH]) and 
(03D for the RI/SMOM scheme leads to zf /SMOM = l/Z^ SMOM . 

The conditions in Eq. fllip for the pseudoscalar and scalar operator can be expressed 
in terms of the bare Green's function and the renormalization constants. The traces over 
the two bare Green's functions become equal in the massless limit in perturbation theory, 
which leads to Zs = Zp. 

In the above discussion the renormalization constants relate the bilinear operators 
renormalized in the RI/SMOM scheme to those in the bare theory which we had assumed 
to satisfy the Ward-Takahashi identities (EJ) and (J7|). Since many lattice formulations 
of QCD break the chiral or flavor symmetries, in general Eqs.flH]) and (|7j) do not hold 
in these (bare) theories. Nevertheless, our renormalization scheme is indeed regulariza- 
tion independent and the Ward-Takahashi identities hold for the RI/SMOM renormalized 
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quantities. The renormalization constants relating the renormalized and bare lattice op- 
erators depend on the regularization of course, so that, for example, Zy and will 
typically be different from 1 in such cases. 



3 Conversion factors: 

Results of the next-to-leading order calculation 

The properties discussed in Section [2] can be used to convert quark masses determined 
through lattice simulations in the RI/SMOM scheme into the MS scheme by computing 
the matching factor C^ SMOM = Z^/Z^ MOM with mf* = C^/ SMOM mg I/SMOM . The 
explicit calculation to determine this conversion factor at one-loop order in perturbative 
QCD will be performed in the next subsections using two different methods, which allows 
us to cross-check our results. 

First, the matching factor cm/smom can ^ e obtained with the help of Eq.Q through 



(c *u/SMOM r i = (C5 /MOM )- 1 - \cf / SMOM lira * Tr 

2 m R^ 12m^ s 



% A A,R 75 



(15) 



ay in 



which will be evaluated in Section 13.11 In analogy to (7 RI / SMOM we define here the con- 
version factor Cf /SMOM = Zp/Z g RI / SMOM for the fermion fields. 

Second, the conversion factor can be related to the renormalization constants of the 
pseudoscalar operator 

7 MS 7 RI/SMOM 1 

.^RI/SMOM _ _ 1 /-i fi \ 

°m ^RI/SMOM 7 ms ~ ^BJ/SMOM v iU J 



and hence 



7 RI/SMOM yMS ~ ^yRI/SMOM 

Jra t> p 

MS _ 1 RI/SMOM _ 1 1 ? n 

m R /-(RI/SMOM m ^ " ^RI/SMOM 7 RI/SMOM mB ' /a **-' V 1 ' > 

In particular in Section I3T21 we will evaluate the conversion factor (j^-/ SMOM ^ w hich con- 
verts the pseudoscalar operator from the RI/SMOM scheme to the MS scheme. The 
matching factor ^Vsmom . g g enera | g aU g e dependent; however, this gauge dependence 
will cancel out with the corresponding gauge dependence in the factor Zp 11 /^ 1 ^ deter- 
mined in the lattice calculation. In the following we will perform the computation in the 
general covariant gauge using the tree level gluon propagator 

.^(V + d-O-gL) (18) 

q 2 + le \ q 2 + it j 

and we will restrict ourselves to the Landau gauge (£ = 0) at the end of the calculation. 
We choose the renormalization scales of both schemes to be equal fx MS = ^ ri /smom^ 
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conversion factors C-m/smom with x e | m? ^ ^ p? y- ^ T | denote always the conversion 
from the RI/SMOM to the MS scheme. 



3.1 The vector and axial- vector operator 

In this section we want to use the vector and axial- vector operator separately to extract 
the matching factor cw/smom f or quark field for the symmetric subtraction point. 
This result is then used in the next step to compute c^/smom ^ e ^ e \p Q f Eq. f |15p . 

Our perturbative computation is performed in dimensional regularization with the 
space-time dimension d = 4 — 2e. For the vector operator case C^ l ' SUOM can be obtained 
by 

^RI/SMOM -j- 1 • '■ 



lim 

m R ^o I2q 2 



Tr 



MS 



(19) 



The calculation of the one-loop QCD corrections to the vector operator, computed between 
two off-shell quark lines, is straightforward and leads to 



^RI/SMOM 



a. 



l-^C F i + 0{a 2 s ). 



(20) 



The symbol Cp denotes the Casimir operator of the SU(3) group in the fundamental 
representation; Cp = 4/3. As expected Eq. fpZO]) agrees with the result in Ref. [6l[7], since 

C RI/SMOM = ^RI'/MOM^ ag ghown in g ection[ 2j 

Similarly one can also derive this result from the axial-vector operator by using 



^RI/SMOM y 



lim 

m a -xi Ylq 1 



Tr 



a u.MS 
% k A,R 75 



(21) 



sym 



For the treatment of 75 in dimensional regularization [2,20] we use a naive anticommuting 
definition of 75 for evaluating the loop integrals, which obeys the equations {75,7^} = 
and t| = 1. This is a self-consistent prescription for the flavor nonsinglet contributions 
considered in this work [21p2] . On the other hand one can use Eq. (jSJ) in order to determine 
^ri/smom an( j then extract ^ I / SMOM anc j ^m/smom f rom Eqs. (|10p . For both we explicitly 



RI/SMOM 



1 as expected. 



confirm that at one-loop order Zy l ^ SMOM _ % 
The conversion factor C^/smom 

can now be computed from the axial- vector operator 
with the help of Eq. (j!5p by determining 



lim 



Tr 



AM.MS 



75 



sym 



Air 



C F (3 + C)C + O(al 



(22) 



with 



C n 



3 V3 



-71 

3 



(23) 
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where is the digamma function ty(x) = T'(x)/T(x) q The matching factor (7^ I '/ MOM 
can be taken from Ref. 0[7] and Cf l/SMOM from Eq-fl^DT). This leads to 

^RI/SMOM = 1 _ ^_ Cf (^ + ^ _ ( g + 1^ + 0{a 2y ^ 

3.2 The pseudoscalar and scalar operator 

In this section we determine the conversion factor (7W/SMOM _ ^w/Smom^_i _ ^ri/smom^-ms 
through the calculation of the pseudoscalar operator. At one-loop order in perturbative 
QCD its computation leads to the decomposition 

m4 „ . ^1,^2] 



with 



Ap, B = A P:B i<y 5 + B P:B i>y 5 —- + C P ^ b — , (25) 



A p ,b = 1 + 7^,1 + • • • , 5ps = -^^,1 + • • • , C P)B = ^-Cpi + (26) 

47T 47T 47T 



where the dots stand for higher order corrections and where we have set pf = q 2 = p 2 ,. 
The quantities b v ,\ and c Pt \ are finite, whereas a Pj i contains 1/e-poles. In the limit of 
massless fermions, considered here, we obtain 6 p l = 0. The matching factor (j^-/ SMOM 
can be obtained from Eq. ffTTj) by evaluating 



^RI/SMOM = ^RI/SMOM ^ J_ Tt 
P 1 mR ^ I2i 



A P,i?75 



(27) 

sym 



The fermion field conversion factor cm/smom . g known and has been discussed in the 
previous Section [XU Since the amplitudes B PB and C P>B in Eq. f[2"5"j) do not contribute 
to the trace in Eq.( lTT|) . this condition depends only on the pseudoscalar amplitude A PtB , 
which is fixed to its lowest order value at the symmetric subtraction point. The explicit 
calculation yields 



^RI/SMOM __ £,RI/SMOm|^ 1 a s 



4+2^+ + (|- 2 -^'Q) +°{^ \ (28 



Inserting Cf' SMOM from Eq.(|20]) and exploiting C^ SMOM = (c^ I/SMOM )- 1 leads to the 
same result as given in Eq. (I24j) . Numerical evaluation in the Landau gauge leads to 

^ri/smom = ^ri/smomx-i = x + ^ .4841391... + O (at) . (29) 

Comparing Eq.® to the RI'/MOM scheme with C^'/ MOM = 1 - f^C F 4 + ... or the 
RI/MOM scheme, (which is in the Landau gauge at one- loop order equal to the RI'/MOM 
scheme), we see that the result in Eq. f[2"9~j) has a smaller one-loop coefficient by almost a 



1 The prime denotes here the derivative. 
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factor of 10. 



In order to study the conversion factor (7W/SMOM f Qr different subtraction points, 
we introduce the parameter oj and fix our renormalization condition for the subtraction 



•'point' 



pi 



—fx and q = —up:. This allows us also to study the limit oj — > 0, 



which results in an exceptional momentum configuration, whereas the limit oj — > 1 gives 
the symmetric one. The result depending on oj is given by 



^yRI/SMOM 



OJ 



+ O (a>) 



(30) 



where the function C (oj) for a; G [0,4] is given by 

'd 4 k 1 



C„(w) 



27r 2 (A; + p±) 2 (k + P2) 2 k 2 



p' 2 =p^ =—^i 2 ,q 2 =—ujfj? 



2i 



Wa/W 




(31) 



and 1^2(2) is the dilogarithm function. In the case oj — 1 one obtains the result of Eq. (j24l) 
with Co(o; = 1) = Co. In order to display the dependence of this result for c^/smom on 
the gauge parameter £, we introduce the one-loop coefficient function c^' RI / SMOM (c<;, £) 
extracted from Eq.fl30]) using the definition: C^/ SMOM = 1 + f^C F c«' RI / SMOM (w, £)• The 
coefficient c^ ,RI//SMOM (a;, £) is plotted as a function of u; in the interval oj G [0,4] for 
different gauges in Fig. Hta) . 




m=0 



co=1 



0.5 1 1.5 2 2.5 3 3.5 4 

CO 

(a) (b) 

Figure 4: (a) shows the one-loop coefficient c C0> ri /smom Q f ^ e matching factor C RI//SMOM 
as a function of oj for different gauges. The value at oj = is the result for the exceptional 
momentum configuration. The nonexceptional configuration is indicated through the 
vertical line at oj — 1. The bold line indicates the Landau gauge (£ = 0), which is usually 
adopted in lattice calculations, (b) shows the exceptional {oj = 0) and nonexceptional 
(oj = 1) configuration as a function of the gauge parameter £. 
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Going from the exceptional (u = 0) to the nonexceptional (u = 1) momentum con- 
figuration leads to a smaller one-loop coefficient in the Landau gauge; however even for 
almost all other gauges the one-loop coefficient becomes smaller as well, except for gauges 
in the small interval £ G ( c , q 4 _ 4 — 3, —3), which is shown in Fig. H^b). The smaller co- 
efficient might indicate that the symmetric configuration is less disposed to infrared effects. 

In analogy to the pseudoscalar operator the computation of the scalar operator leads to 
the conversion factor (j^' SMOM by employing the renormalization condition in Eq. fllip . As 
expected the one-loop result for the matching factor C^ I//SMOM = ((^VSmom-j-^ ex ^ rac ^ ec | 
from the scalar operator, is equal to the result obtained from the pseudoscalar one. 



3.3 The tensor operator 



The matching factor converting the Green's function of the tensor operator from the 
RI/SMOM to the MS scheme can be obtained from Eq. ffT2l in complete analogy to the 
other operators. It is given by 



^RI/SMOM 

Lsnr 



'T 



:i-o 



Co 
2 



(32) 



u t, j c j ^yRI/SMOM 

where we have defined C T 
Landau gauge leads to 



7 MS /7 RI/SMOM 

Zlrr 



a 



RI/SMOM 



The numerical evaluation in the 



1 -C F 0.1613797. 

4:71 



(33) 



For the RI'/MOM scheme the tensor operator has been evaluated up to three-loop order 
in Ref. [7]. The conversion factor at one-loop order is found to be proportional to the 
gauge parameters £. This contribution is therefore zero in the Landau gauge. 

All conversion factors discussed in Section [3] are summarized in Table [2J The matching 
factors for the scalar and pseudoscalar operator are equal to the inverse of the mass 



conversion factor i/cM/smom = c 



RI/SMOM 
S 



C 



RI/SMOM 



£<RI/SMOM _ 


1 - f^CV 0.4841391.. 


+ o (a?) 


£<RI/9MOM _ 


1 + O (a 2 s ) 




^yRI/HMOM 


1 - f^C F 0.1613797.. 





Table 2: Summary of the matching factors for the mass and fermion field conversion 

^RI/SMOM (■}{] SMON ' 

in the Landau gauge. 



^ri/smom ) c ri/smom^ as well as for the conver sion of the tensor operator ((7? /SMOM ) 
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4 Summary and Conclusion 



We provide the framework and concepts for renormalizing the quark bilinear operators in 
a MOM scheme (RI/SMOM) with a symmetric subtraction point which has no channels 
with exceptional momenta. This generally suppresses the infrared chiral symmetry break- 
ing effects compared to the standard RI/MOM (or RI'/MOM) scheme in which there is an 
exceptional channel (with zero momentum). An exception is the vector current for which 
the RI/MOM scheme satisfies the Ward-Takahashi identity also at low values of p 2 . We 
demonstrate that the chiral Ward-Takahashi identites(for degenerate masses) are satisfied 
nonperturbatively, and thus Zy = 1 = for all values of p 2 , in the RI/SMOM scheme. 
We calculate the matching factors relating operators renormalized in this scheme and the 
MS schemes at one-loop order in perturbation theory. The one-loop coefficients are given 
in Tableland we note that they are small. In particular, for the quark mass the coefficient 
is much smaller than that between the RI/MOM(RI'/MOM) and MS schemes which, if 
confirmed at higher orders, would lead to a significant reduction in the uncertainty on the 
calculated value of the quark mass. 

Nonperturbative renormalization of operators in lattice QCD using the RI/MOM (or 
RI'/MOM)scheme has been successfully implemented for many years. The evaluation of 
matrix elements in the RI/SMOM renormalization scheme in lattice simulations is equally 
practicable and in view of the advantages explained above we strongly advocate its use. 
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A Alternative projectors for the vector and 
axial-vector operator Green's functions 

In general one can also use other projectors than those of the RI/SMOM scheme as 
defined in Section [2] in order to define a scheme with a symmetric subtraction point. 
The general structure before taking the trace with projectors of the one-loop corrected 
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amputated Green's functions of the operators in the MS scheme for massless quarks with 
the momenta p\ = p\ = q 2 are given by 



As 
Ap 



A s l + C s — — dij, 



A P i-f 5 + C P i>y 5 



AM), 



Ai^A+Ai^A , ^ Ai 



V 



— + L> V 2 



D f 



Q 



where the indices % and j denote color indices and the coefficient functions read 



(34) 
(35) 
(36) 

IM37) 
(38) 



A s ,p 
Cs,p 
A v ,a 

By,A 
Cv,A 
D V ,A 

A T 
Ct 



i + ^c f 

An 



4 + 3 log 



-? 2 



-Q 



\_Co_ 
2 j 2 



\_Co 

2 ] 3 



-^f [(1-0(^0-2) -2], 
-g^[Co + 2 + e(C -l)], 
£f [(I" O(^o-l) -1], 



1-^(1-0 [3+ log 

-^pcb-(i-0], 



-9 



2 j 



4tt 3 



1-0 [2-Co]. 



(39) 
(40) 
(41) 

(42) 
(43) 
(44) 
(45) 

(46) 
(47) 



An example of a second possible choice for the projectors is the use of the projectors of 
the RI/MOM scheme for the amputated Green's function of the vector and axial- vector 
operator 



lim — Tr 

m R ^0 48 



A£ )fl (pi,p 2 )7/i 



1, 



sym 



lim — Tr 

m R ^0 48 



A A ii? (Pl,P2)757M 



1, (48) 



sym 
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in the renormalization conditions with a symmetric subtraction point together with the 
conditions of Eqs.( JTT]) . One also has to modify the conditions of Eqs.(@| 



lim — < Tr 

m R ->0 48 



Tr 



7 



lim 



s r \p) 



+ Tr 



%1 



dq a A V,R 



-1, 



-c 2 2 



sym _ 



sym 



(49) 



(50) 



to maintain the Ward-Takahashi identities of Eqs. ([6]) and (J7J) for renormalized quantities. 
This leads to a wave function renormalization factor Z q which is different from the one 
of the RI/MOM or RI'/MOM scheme. For this reason the projectors used in Eqs.lflOl of 
Section [2] have the advantage to produce the same well-known renormalization constant Z q 
like in the RI'/MOM scheme. With the conditions of Eqs. flTT]) and (j4"8"]) - (l5T)]) one obtains 
in this RI/SM0M 7 scheme the following conversion factors: 



RI/SMOM^ 



4tt 



2 I 3 V3 



-71 

3 



+ a 



RI/SMOM^ 

Up 



5+ r 2 



04 



+0 (of) . 



(51) 



(52) 



The numerical evaluation of the resulting mass conversion factor in the Landau gauge 
reads 

ri/smoivt a. 
C m " = 1 - ^Cjrl.4841391 . . . . 

4:71 



(53) 



B Anomalous dimensions 



In order to evaluate the mass in the RI/SMOM scheme at different scales the mass anoma- 
lous dimension y*&/ SMOM [ s required. It is defined by 



7r, 



dlogm(fi) 
d\og (/i 2 ) 



Im Im 



71 



71 



0{al). 



The result up to order a 2 s in the RI/SMOM scheme reads in the Landau gauge 



(0),RI/SMOM = (0),MS 
Im Im ) 



(i),ms _ P^_ n ji), 



with C W' RI / SMOM (1,0) = -0.4841391... and the j3 function defined through 



(54) 



7 (l),RI/SMOM = 7 (1),MS _ ^_ CfC (1),RI/SMOM (1; ^ ^ 



(56) 



The first expansion coefficients for the mass anomalous dimension in the MS scheme and 
the MS (3 function are given by 



(O)-MS = | Cfj 



=\(fC A -lT F n f ) 



~(1),MS 
Im 16 
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where Ca is the Casimir operator in the adjoint representation of SU(3) and n/ is the 
number of active fermions. The symbol Tp denotes the normalization of the trace of the 
SU(3) generators in the fundamental representation, conventionally chosen as 1/2. 
For the RI/SMOM 7 scheme, defined in Appendix El the two-loop mass anomalous di- 
mension is given by 

(1),RI/SMOM 7m _ — 0(0) (l),RI/ S MOM 7M 
7m — 7m —^Ftm , {O I ) 

(l),RI/SMOM, 

with Cm = -1.4841391 ... as given in Eq. fl53j). 

Similarly the anomalous dimension ^m/smom _ 2_^i2S-^. can De defined, which is equal 
to the one in the RI'/MOM scheme, 7 f/ SMOM = ^bi'/mom^ and can be found in Refs 
up to order a^. For completeness we give here the result up to 0(a 2 s ) which in the Landau 
gauge is the same as in the MS scheme 



7 ki/smom = ^ _ ^ CfCa + l _c F T F n f ) + 0{al). (58) 



In the RI/SMOM scheme, defined in Appendix |A] the order a 2 coefficient of the anoma- 

RI/SMOM 

lous dimension ^ q M in the Landau gauge reads 

We define the anomalous dimension 7^ of the tensor operator by 

RI/SMOM d\ogZ T (0) { a A (1) ( ®s\ 2 . , nf 3 n , Rn s 

In the RI/SMOM scheme it is given in the Landau gauge by 

7 (0),RI/SMOM = 7 ( ),MS ) 7 (l),RI/SMOM = ^D.MS _ ^^D.RI/SMOM (Q) ? 

where we have introduced the one-loop coefficient function Cj^' RI ^ SMOM (£) extracted from 
Eq.(E2D using the definition C? /SMOM = 1 + ^C F c^ m/SM0M (0 , with c?' R1/SMOM (0) = 
—0.1613797.... The MS anomalous dimension 7 ^ s is known from Refs. [7|l2~3|i24"] and reads 

(o),ms l^y (i),ms 1 / 19 2 257 26 \ 

7t = ^<^f, 7t = ^ (."Y^ + T8~ CfCa ~ "9 ' C * FU} ) ■ ( 6 ^ 

Since the renormalization constants of the pseudoscalar and scalar operator are related 
to the mass renormalization constant, the anomalous dimensions of the pseudoscalar and 
scalar operator follow from the mass anomalous dimension in Eq.( l55l) . 
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